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Abstract. The lattice vibrations are determined in the quasi-harmonic approximation for elemental Au and Cu to twice their 
normal density by first-principles electronic band-structure calculations. It is found for these materials that the important 
moments of the phonon density of states can be obtained to high accuracy from short-ranged force constant models. We 
discuss the implications for the Gruneisen parameters on the basis of calculated phonon moments and their approximations 
by using bulk moduli and Debye- Waller factors. 
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Accurate and reliable high-pressure standards are fun- 
damental for the study of matter under extreme condi- 
tions and in the earth sciences. With the advent of neu- 
tron and x-ray diffraction measurements up to several 
hundred GPa in diamond anvil cells it becomes impor- 
tant to be able to determine precisely the pressure. At 
high pressures the equation of state (EOS) of elemental 
metals are commonly used as pressure scales. However, 
these scales are based on either extrapolation of low- 
pressure data or reduction from the Hugoniot to room 
temperature. The reduction of the Hugoniot curve onto 
P — V isotherms requires model assumptions about the 
Gruneisen parameter, which are experimentally not well 
founded. Our ability to calculate accurate moments of 
the phonon distribution from first-principles electronic 
structures enables us to strongly constrain the Gruneisen 
parameter and to establish an accurate EOS and high- 
pressure standard for elemental gold and copper. This 
is an important step toward high-precision high-pressure 
experiments. 

Since recent experiments raised concerns about the 
consistency of these standards fT], we revisited the prob- 
lem of high-pressure standards using a semi-empirical 
approach combined with first-principles electronic struc- 
ture calculations |2]. This method is different from the 
tight-binding approach that was earlier applied to sohd 
copper |,3:]. We determined the lattice vibrational con- 
tribution to the EOS in the quasi-harmonic approxi- 
mation of elemental face centered cubic Au and Cu 
by first-principles electronic structure calculations. We 
found that the important moments of the phonon density 
of states can be obtained to high accuracy from short- 
ranged force constant models. 



In this paper, we present the implications for the 
Gruneisen parameter based on those calculations and 
their approximations by using bulk moduli and Debye- 
Waller factors which, in principle, can be obtained 
from ultrasound and neutron/x-ray diffraction measure- 
ments. The accuracy of our results is based on the ac- 
curacy of first-principles electronic structure calcula- 
tions, which were performed in the local-density approx- 
imation (LDA) and generalized-gradient approximation 
(GGA) of density functional theory. From these results 
elastic moduli and zone boundary phonon frequencies 
were obtained. Fitting simultaneously a short-ranged, 
second nearest-neighbor Born-von Karman (BvK) force 
matrix model to the zone boundary frequencies and the 
elastic moduli enabled us to generate phonon disper- 
sion curves in the entire Brillouin zone and to compute 
any phonon moment needed, which is typically accurate 
within a few percent. 

At high temperatures, i.e., in the classical limit, the 
temperature dependence of the pressure of a solid is 
dominated by the contribution from lattice vibrations, 
P^*'. Its temperature derivative at constant volume, 
{dP^'^/dT)v, is proportional to the Gruneisen parame- 
ter 7, which is defined by 
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where the phonon frequencies o\ are functions of vol- 
ume only and the summation is over all eigenmodes. 
Eq. Q defines the logarithmic phonon moment (Oq. 

For practical reasons one often uses an interpolating 
Debye phonon model for calculating the EOS or for ana- 



TABLE 1. Bulk modulus B, log-moment (Oq, thermal mean-square displacement {u ), and 
Griineisen parameter 7 for Au and Cu at ambient conditions were obtained from the theoretical 
data in Figs. 1(a) and 2(a) using Eq. 0. Experimental results are listed where available. 
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lyzing diffraction data, instead of the more elaborate lat- 
tice dynamical models. Since in a Debye phonon model 
the Debye frequency Ofe is identical to all other phonon 
moments, the approximation (Od ~ 0„, with n > —3, is 
widely used for computing lattice vibrational properties 
and Griineisen parameters. The moments a)_2 ™d a)_3 
are related to the Debye- Waller factor at high tempera- 
tures and the sound velocity, respectively |7]. Using the 
theoretical relationships for the thermal mean-square dis- 
placement of an atom, (m^) 0= T jaP^^, and the sound 

speed, c °^ o_3 °^ vBV^, we find the following ap- 
proximations to the Mie-Griineisen theory for the high- 
temperature Griineisen parameter. 
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The derivation of Eq. Q requires a constant Poisson 
ratio (Slater approximation). Moruzzi et al. studied 
extensively its application to the 4d transition elements 
and found good agreement with experiment at ambient 
conditions. The advantage of using Eq. (|3 for estimat- 
ing the Griineisen parameter is that it is readily accessi- 
ble from ultrasound and diffraction measurements. On 
the other hand, the expression for the thermal param- 
eter in Eq. has not seen much application, because 
of the extreme difficulties of measuring accurate thermal 
mean-square displacements in high-pressure diffraction 
experiments |10|. Since these different approximations 
of the Griineisen parameter emphasize different phonon 
frequencies in the Brillouin zone compared to the log- 
moment cOq, we do expect to find deviations from 7 
by at least several percent, when using these approxi- 
mate formulas, reflecting the differences between differ- 
ent phonon moments. 

The calculation of the phonon moments co„ requires 
knowledge of the phonon frequencies for all k points in 
the Brillouin zone |7J. The direct first-principles calcu- 
lation of frequencies on a dense k mesh is computation- 
ally intensive, while the phonon dispersions of most ele- 
ments can be easily parameterized using lattice dynami- 
cal models like a generalized BvK force matrix model. 



Often it suffices to use a short-ranged force model to 
compute the low order moments accurately within a few 
percent 1 11 ] . In particular, for elemental Au and Cu the 
log-moment coq is converged to less than 1 % with a 2nd 
nearest-neighbor interatomic shell model at ambient con- 
ditions. Thus, we chose to calculate four zone boundary 
phonon frequencies coiTesponding to the transverse and 
longitudinal eigenmodes at the X and L points of the Bril- 
louin zone. These are computed with standard frozen- 
phonon methods. Additionally, three elastic moduli are 
computed using the method by Soderlind et al. lll2ll . We 
fitted these results to a 2nd nearest-neighbor BvK force 
model, which then allows the evaluation of the frequen- 
cies (Ok in the entire Brillouin zone. Details of this calcu- 
lation will be published elsewhere f?|. 

We fitted the theoretical (Oq, B, and (m^) results to a 
functional form that gives a realistic V-dependence of the 
Griineisen parameter, which has been used in many EOS 
calculations [13], 



7(y) = f +Ai (V/Vo) +A2 {V/Vof 
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where Vq is the volume at ambient pressure, 7°° is the 
infinite density limit of 7, and Aj and A2 are fit parame- 
ters. Recently, it has been argued for the value 7°° — 1/2, 
instead of the commonly used 7°° = 2/3 1 14]. However, 
our results are insensitive to this difference and hence 
we chose 7°° = 2/3. The corresponding fitted values are 
listed in Tabled 

In order to test the robustness of the calculated 7 values 
using Eq. (0}, we also fitted too, B, and (m^) to the widely 
used expression 
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where the fit parameters y{Vq) and q are listed in Table 
121 see also the inserts of Figs. 1(b) and 2(b). 

In Figs. (1) and (2) we show the normalized log- 
moments, bulk moduli, and thermal mean-square dis- 
placement parameters of elemental Au and Cu based on 
our first-principles calculations. To emphasize the very 
similar scaling behavior of these lattice dynamical prop- 
erties, we normalized their values by their corresponding 
ambient condition values. 

We used the functional forms for the Griineisen pa- 
rameter given in Eqs. Q and (|5jl to integrate Eqs. 



TABLE 2. Fit parameters and coq, B, and (k^) for Au and Cu at ambient conditions were obtained from fitting the 
theoretical data in Figs. 1(a) and 2(a) using Eq. jsj- The reference values y^ef ars from 1 7] and were derived using the 
thermodynamic definition of y= {V /Cv)idP^''^ /dT)v, with the specific heat Cy- 
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FIGURE 1. (a) Normalized (to ambient conditions) log- 
moment, bulk modulus, and thermal mean-square displace- 
ment of Au at r = 296 K from electronic structure calcula- 
tions at V/Vo = 0.5,0.6,0.7,0.8,0.9, 1.0, 1.1. (b) Correspond- 
ing Griineisen parameters using Eq. y(Vo) = 2.95 Q] 
(cross) is shown for reference. Insert: Corresponding Griineisen 
parameters using Eq. Jsj- 




FIGURE 2. (a) Normalized (to ambient conditions) log- 
moment, bulk modulus, and thermal mean-square displace- 
ment of Cu at r = 296 K from electronic structure calcula- 
tions at V/Vo = 0.5,0.6,0.7,0.8,0.9,1.0,1.1. (b) Correspond- 
ing Griineisen parameters using Eq. y(Vo) = 2.02 1"^ 
(cross) is shown for reference. Insert: Corresponding Griineisen 
parameters using Eq. Js}- 



(O through 0. This allowed us to fit the theoretical 
log-moments, bulk moduli, and thermal parameters and 
to extract the fitting parameters necessary for calculat- 
ing the corresponding Griineisen parameters shown in 
Figs, (lb) and (2b). 

We found for the volume compressions considered 
here that 7 of Au and Cu are quite well described by 
either expression @ or (|5}. The differences in the cal- 
culated 7 values are only a few percent. A detailed dis- 
cussion of the accuracy of the different functional forms 
for 7(y) and its consequences for the equation of state 
and the Hugoniot is given in Ref. |2,]. In the case of Au, 
the approximate expressions for the Griineisen parameter 
in Eqs. (|2j and Q give good agreement with the correct 
high-temperature 7 obtained from Eq. ([O- The deviations 
of Yb and yow from 7 are mostly less than 8%. In the case 
of Cu, the deviations of Yb and Jdw from 7 are generally 
bigger than for Au, but never more than 15% for com- 
pressions in the range 0.5 < V/Vo < 1-0. A possible ex- 



planation of the larger deviations for Cu may be that the 
Debye temperature and equivalently the log-moment of 
Cu are almost twice as large as for Au. Since the typical 
temperature of lattice vibrations of Cu, hcoo / ks ^ 330 K, 
is slightly above room temperature (300 K), one may 
have to go beyond the classical approximation for ac- 
curately calculating the Griineisen parameter at ambient 
conditions. 

Comparing jow and Jb with 7, we find in the range of 
compressions 0.8 < V/Vo < 1-0 that the Debye-Waller 
approximation yow results generally in better agreement 
with 7 than the approximation using the bulk modulus 
Yb- This situation is reversed below V/Vo ~ 0.8. At such 
high compression 7^ is in very good agreement with 7 
of Au and in good agreement for Cu, while 701^^ devi- 
ates the most. A simple explanation of this very different 
behavior of Yb and Ydw at high compression, i.e., below 
V/Vo ^ 0.8, is due to the drastic stiffening of the phonon 
frequencies and the simultaneous increase of the bulk 



modulus with decreasing V. At such high pressures and 
room temperature the excited phonons probe mostly the 
linear part of the phonon dispersion. The slope of the dis- 
persion near the zone center is crudely proportional to the 
bulk modulus, while the Debye- Waller factor averages 
all frequencies weighted by the temperature dependent 
occupation factor of each mode. Therefore, we expect 
Eq. (|3} to fail when the temperature becomes comparable 
to the high-temperature Debye- Waller phonon moment. 
Roughly speaking, for T of the order of the log-moment, 
ksT ^ TicOq. 

In summary, we have successfully computed with high 
accuracy phonon moments of elemental gold and cop- 
per from first-principles electronic structure calculations 
combined with a Born-von Karman force model. From 
the logarithmic phonon moments we calculated the vol- 
ume dependence of the Griineisen parameter up to twice 
the normal density of Au and Cu at ambient conditions. 
Comparing the Griineisen parameter with approxima- 
tions based on the bulk modulus and the Debye- Waller 
factor, we found that for low compression, 0.8 < V /Vb < 
1.0, the approximation using the thermal mean-square 
displacement is more accurate than the one using the bulk 
modulus, while at higher compression the bulk modu- 
lus gives generally better agreement with 7. Therefore, 
a combination of Jb and Jdw^ which can be obtained 
from ultrasound and diffraction measurements, can give 
an estimate of the volume dependence of the Griineisen 
parameter 7 within approximately 10%. This provides a 
useful alternative for determining the Griineisen parame- 
ter besides using the thermodynamic relation, which de- 
pends on the knowledge of the bulk modulus, thermal 
expansion, and specific heat. 
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